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Abstract 

We consider a family of integro-differential equations depending upon a parameter b as 
well as a symmetric integral kernel g{x). When 6 = 2 and g is the peakon kernel (i.e. 
g{x) = exp(— |x|) up to rescaling) the dispersionless Camassa-Holm equation results, 
while the Degasperis-Procesi equation is obtained from the peakon kernel with 6 = 3. 
Although these two cases are integrable, generically the corresponding integro-PDE is 
non-integrable. However, for 6 = 2 the family restricts to the pulson family of Fringer 
& Holm, which is Hamiltonian and numerically displays elastic scattering of pulses. On 
the other hand, for arbitrary b it is still possible to construct a nonlocal Hamiltonian 
structure provided that g is the peakon kernel or one of its degenerations: we present 
a proof of this fact using an associated functional equation for the skew-symmetric 
antiderivative of g. The nonlocal bracket reduces to a non-canonical Poisson bracket 
for the peakon dynamical system, for any value of 6 7^ 1. 



1 Introduction 

Here we will consider a class of integro-partial differential equations (IPDEs) of the form 

mt + unix + bUxTn = 0, u = g * m, (1.1) 

where g* denotes convolution with a symmetric integral kernel g{x) defined on the real line 
R: 

/oo 
g{x — y)m{y) dy. 
-00 

There are some distinguished special cases of this equation, the first being the dispersionless 
version of the integrable Camassa-Holm equation, 

— Uxxt + 3nn^. = 2uxUxx + uu^xx^ (1-2) 
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which is the b = 2 case of when the kernel is chosen to be the Green's function for the 
Helmholtz operator on the hne: 

= ^e"l''l, m = u-u^^. (1.3) 

With the inclusion of additional linear dispersion terms the equation p.2j) was derived as an 
approximation to the incompressible Euler equations, and found to be completely integrable 
with a Lax pair and associated bi-Hamiltonian structure H] (for another derivation, see 
|31j). In fact the Camassa-Holm equation fits into the framework of hereditary symmetries 
and recursion operators described earlier by Fokas and Fuchssteiner [22] • A remarkable 
discovery of jHl EI was that in the dispersionless limit ()1.2j) . the solitons of the Camassa- 
Holm equation are peakons, given by a superposition of an arbitrary number of peaks, 

1 ^ 

u{x,t) = -J2pAt)e~^'-''^'\ (1.4) 
i=i 

(See also [21 CHI CI]-) The A^-peakon solutions (jl.4p are weak solutions with discontinuous 
first derivatives at the positions qj of the peaks; qj, pj are canonical coordinates and momenta 
in an integrable finite-dimensional Hamiltonian system; the interpretation of weak solutions 
is discussed further in ^2). The A^-peakon interaction was obtained in jHj, with the explicit 
peakon-antipeakon interaction being given in |3]. The stability of peakons in the Camassa- 
Holm equation was proved in ^Sl ■ 

In a recent application of the method of asymptotic integrability to a many-parameter 
family of third order equations with dispersion, Degasperis and Procesi J7] found that only 
three equations passed the test up to third order in the asymptotic expansion, namely KdV, 
Camassa-Holm and one new equation. After rescaling and applying a Galilean transforma- 
tion, the new equation may be written in dispersionless form as 

Ut - U,j;^t + 4:UU^ = Su^Uxx + UUxxx- (1-5) 

Henceforth we refer to p.Sj) as the Degasperis-Procesi equation, and observe that it is the 
b = 3 case of the equation (jl.lll with g{x) = e~'^'/2. In a recent article ^H] we proved the 
integrability of the Degasperis-Procesi equation by constructing a Lax pair, and derived two 
infinite sequences of conservation laws. 



u(x,t) u(x,t) 




Figure 1: The elastic scattering of two peakons in the integrable Degasperis-Procesi equation 
on a periodic domain. This corresponds to b = 3 and g(x) = e~'^'/2 in 

In order to better understand the common properties of the two integrable equations 
(jl.2j) and (jl.Sp . we have found it convenient to consider the following one-parameter family 



of partial differential equations (PDEs): 



Ut - Uxxt + {b + l)uu^ = bur^u^^ + uu^,,^. (1.6) 

Since it arises from p.l|) when the peakon kernel ()1.3|) is chosen, we refer to fll.6|) as the 
peakon 6-family of PDEs. Remarkably, every member of this family admits multi-peakon 
solutions ()1.4|1 . with a non-canonical Hamiltonian structure for the dynamics of qj, pj found 
in PHI- The two-body dynamics is integrable for any b (see Figure^, and for all values of 
b > 1 the peakons appear to be numerically stable and dominate the initial value problem 
|24j . Furthermore, with the addition of linear dispersion {u^ and u^xx terms) it has been 
shown that not only the Camassa-Holm equation 120] but also the whole 6-family ()1.6|1 (apart 
from b = —1) belong to a class of asymptotically equivalent shallow water wave equations 

In we presented a Lagrangian formulation for the 6-family ()1.6|) . with a Legendre 
transformation leading to the Hamiltonian structure 

where for 6 7^ 1 

B = -{bmd, + m,){d,-dl)-\bmd^ + {b-l)m,), H = J mdx. (1.8) 

(In the case 6 = 1 it is necessary to take H = J mlogmdx instead.) The skew-symmetric 
operator B in p.8|) was first obtained in ^Hl for the case 6 = 3, as a second Hamiltonian 
structure for the Degasperis-Procesi equation ()1.5|) . However, a proof of the Jacobi identity 
for this operator was lacking until the work of one of us with Wang where the trivector 
formalism of Olver jHSj was applied to this problem. The first two Hamiltonian structures 
of (1131) are 

Bo = -8,(1 - dDiA - dl), B, = 5|5=3 = -Qm"^d^m^/\d^ - dlY^m^'^'d^m'/^ (1.9) 

while for the Camassa-Holm equation (jl.2|l the bi-Hamiltonian structure arising from the 
Lax pair in [HI d is 

Bo = -d,{l - dl), 5i = -{md^ + d.,m), 
and applying the recursion operator R = BiB^^^ to Bi gives the nonlocal operator 

B2 = BiBq ^Bi = B\b=2- 



Moreover, it is proved in j29j that 6 = 2,3 are the only parameter values for which the 
operator B given by ()1.8p is compatible with another operator with constant coefficients. 
The integrable cases 6 = 2,3 can also be isolated by the perturbative symmetry approach 
|33] . by the Wahlquist-Estabrook prolongation algebra method [SH], or by using a reciprocal 
transformation and then applying Painleve analysis [HH] . 

In the special case of = 2 with an arbitrary symmetric kernel g{x), the equation ()1.1|1 
has the Lie-Poisson structure 

5H 

rrit = {m, H}lp := -{md^ + d^m)— = -adlm, (1-10) 

dm 



with the quadratic Hamiltonian 



H = - / mq * m ax, - — = u. 
2 J 5m 

Equivalently the equation p.lOj) corresponds to the Euler-Poincare equation for geodesic 
motion on the diffeomorphism group of the real hne, assuming that g (the co-metric) is the 
Green's function for a positive definite operator (see [23 for more details, and see [211 for 
an extension to higher dimensions). In the periodic case the Camassa-Holm equation is the 
equation for geodesies on the diffeomorphism group of the circle jHl], leading to a proof of 
the corresponding least action principle • 

In the case 6 = 2 given by ()1.10|) the equation admits special solutions in the form of a 
finite superposition of pulsons whose shape is determined by g, that is 

N N 

u = '^Pj{t)g{x-qj{t)), m = ^pj{t)6{x-qj{t)). (1.11) 
3=1 j=i 

The finite-dimensional dynamical system for the motion of pulsons on the line is a canoni- 
cal Hamiltonian system for geodesic motion on an A^-dimensional manifold with coordinates 
qj, j = 1, . . . , N and co- metric g{qj — qk), derived from the Hamiltonian 

h=^^PjPk9{qj - Qk)- (1-12) 

The canonical Poisson bracket 

{qj,Pk} = Sjk {qj,qk} = = {pj,Pk} (i.i3) 

can be obtained by substituting the ansatz (jl.llj) for m into the Lie-Poisson bracket 

{m{x),m{y)}Lp 

defined by the Hamiltonian operator —{mdx + dxm), and similarly h in p.l2|) is derived 
from H in p.lOj) by substitution. Furthermore, although the equation ()1.10|) is generically 
non-integrable, numerical studies in [221 provide strong evidence that the pulson solutions 
(jl.lip are stable and are produced by arbitrary smooth initial data, just as for the integrable 
case of peakons in the Camassa-Holm equation. We have proved non-integrability for a 
particular fifth order equation in this class [23], although numerically the pulsons appear to 
scatter elastically (see figure 2). 

One might wonder why we have chosen to discuss this subject in a birthday volume in 
honour of Francesco Calogero. In fact Calogero has already studied the problem of finding 
Lax pairs for integrable finite-dimensional systems with geodesic Hamiltonians of the form 
()1.12|) . and obtained functional equations in terms of the function g{x) and the components 
of the Lax matrix (see jH] and references). Here we are concerned with a slightly different 
problem. It turns out that the pulson ansatz p. Ill) works equally well for the whole family 
of IPDEs (II. Ij) . and leads to the dynamical system 

dq- ^ dp- ^ 

-j^ = ^Pk9{qj - qk), -j^ = -{b-l)^PjPk9'{qj -qk)- (1-14) 

k=l k=l 

Clearly this takes the canonical Hamiltonian form, with the Hamiltonian h given by ()1.12p . 
only for the case b = 2. However, in the particular case 6 = 3 and g{x) = e~'^''/2 that 
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Figure 2: Pulson solutions of equation ()1.10|) with g{x) = 2e~'^'' — e"^'"^'. emerge from 

a Gaussian of unit area and width cr = 5 centered about a; = 33 on a periodic domain of 
length L = 100. The fastest pulson crosses the domain four times and collides elastically 
with the slower ones. 

corresponds to the integrable Degasperis-Procesi equation, we have derived a matrix Lax 
pair for the system ()1.14|) . and therefore in this case the system must be integrable with 
respect to some Poisson structure. Requiring the existence of a certain type of Poisson 
structure for the IPDEs (jl.lj) leads us to the functional equation 

G'ia)iG{(3) + G(7)) + G"(/3)(G'(7) + G'(a)) + G'(7)(G(a) + Gi(3)) = (1.15) 

when a + /? + 7 = 0, for the antisymmetric kernel 

G{x) = [ g{t)dt. (1.16) 
Jo 

Apart from some limiting cases, the general odd solution of (jl.l5|) turns out to be 

G{x) = A sgn(x) (^1 - e"-^!""! j , whence g{x) = AEe"^'^' {A, B constant). 

(The derivation of the general solution of (|1.15|) is presented in an Appendix provided by 
H.W. Braden and J.G. Byatt-Smith.) Thus the peakon 6-family ()1.6|) is picked out of the class 
(jl.ll) by the functional equation p.lSj) . Moreover, for any h the corresponding Hamiltonian 
operator, given by p.8|) . reduces to a Poisson bracket for the iV-peakon dynamical system, 
first presented without proof in PU]- In the next section we derive the functional equation 
from the Jacobi identity for the operator (jl.8p . and then proceed to show how this leads to 
a suitable bracket for the peakons. 

Historical Note. This paper is dedicated to our friend, Francesco Calogero. This work 
parallels some of his earlier considerations on geodesic flows and the determination of in- 
tegrability for dynamical systems by using functional equations. Calogero's fundamental 
works have been and remain incredibly influential and important to the field of integrable 
systems. Francesco himself is still developing the mathematics associated with this area. 
(Just look over the 200+ papers on MathSciNet that he has written!) His contributions will 
surely influence the next generation of people working in this RrGci, cis they have with the 
current generation. Of Francesco's work, the references that are most directly relevant to 
the topics discussed in the present paper are jHl El • 



2 Derivation of the functional equation 

Let us start from any equation (jl.lj) in the pulson 6-family. The key is to introduce, for any 
kernel g, the operator 

B = -{bmd, + m,)G{bmd, + {b - l)m,), (2.1) 
where the operator G acts according to 

/oo 
Gix-y)f{y)dy, 
-oo 



with G being the antisymmetric kernel ()1.16jl . Thus B is skew-symmetric and for b 1 
satisfies 

S H 1 /'OO 

mt = B— , H = — / mdx. 

6m {b-l)J_^ 

Now we can calculate the bracket defined by the operator B, namely 

{F,G) = / —B—dx, 
J_oo dm dm 

which gives 

{m{x),m{y)} = {G{x - y)m^{x)m^{y) + bG'{x - y){m{x)m^{y) - m{y)m^{x)) 



-b^G"{x -y)m{x)m{y)y (2.2) 



Although (j2.2|) is antisymmetric (since G is odd), it only defines a Poisson bracket if the 
Jacobi identity is satisfied, i.e. 

{{m(x), m(?/)}, m{z)} + cyclic = 0. (2.3) 

To calculate ()2.3p . we evaluate 

{m{x),m{y)} = (g{x - y)[6'{x - s)m^{y) + 6'{y - s)m^{x)] 



Sm{s) 

+bG\x — y)[S{x — s)mx{y) — S{y — s)m^(x) + 6'{y — s)m{x) — 6\x — s)m{y)] 
—b'^G"{x — y)[5{x — s)m{y) + 6{y — s)m{x)] 
and then substitute into 

{{m{x),m{y)},m{z)} = / - — —{m{x),m{y)}B{s)6{z — s)ds 

J_oo dm{s) 

(where B{s) indicates the operator acting in the independent variable s). 

In the Jacobi identity ()2.3p many terms cancel to leave only terms cubic in m and m^., 
whose coefficients must each vanish if the identity holds for arbitrary m. This leads to 
functional equations for the odd function G, but remarkably (apart from overall prefactors) 
these coefficients are independent of the parameter b. Thus, setting 

a = X — y, P = y — z, 7 = 2 — x, whence a + /5 + 7 = 0, 



we find the coefficient of mx{x)mx{y)mx[z) yields tlie functional equation p.l5|) . Tlie coef- 
ficient of m{x)mx{y)mx{z) gives 

G\l3){G\a) - G'in)) - G"(7)(G(a) + + G"(«)(G(/3) + ^(7)) = 0, (2.4) 

and there are two other identities like ()2.4j) obtained by cyclic permutations of x, ?/, z which 
induce cyclic permutations of a,/?, 7 (with the three identities summing to zero). Similarly 
the coefficient of mx{x)m{y)m{z) yields 

G"\(3){G{a) + G(7)) + G"{-i){G\a) - G'm - G"{a){G\(3) - G\^)) = 0, (2.5) 

and there are two other identities obtained by cyclic permutation of a, /9, 7 in ()2.5p . Finally 
the coefficient of m{x)m{y)m{z) requires 

G"\a){G\(3) - G"(7)) + G'"(/3)(G'(7) - G\a)) + G"\^){G'{a) - G'm = 0. (2.6) 

There are several things to observe about these functional equations. Firstly, ()2.4|) and 
()2.5p are direct consequences of ()1.15p . obtained by differentiating. The final functional 
equation ()2.(jj) is also a consequence of the rest. Hence the Jacobi identity is satisfied by the 
operator B if and only if the antisymmetric kernel G satisfies the functional equation 
We can prove the following: 

Proposition. Suppose that G{x) is an odd solution of the functional equation il.l^) with 
continuous second derivative G"{x) for all x G R\{0}, and suppose lim^^Qj-G" (x) is finite. 

Then either G{x) = Ax or G{x) = Asgn{x) ^1 — e~^l^'l j for some constants A, B. 

Sketch of proof. If G is odd and continous second differentiable then from the limit a — 
in ()1.15|) we see that G(0) = 0, while differentiating (jl.lSp with respect to a yields fl2.4j) . 
Taking the limit a — >■ in (|2.4p with (3 fixed gives the differential equation 

G"{(3)G{P) - G\(3f + G\Q)G\f3) = 0, 

which integrates to 

G'(/3) = G'(0) + i^±G(/?), 

where the integration constants K± can be different for /? > and /? < 0. This is easily 
integrated to give the solution for G. 

Remarks. A more detailed derivation of the general solution (not necessarily odd), with 
weaker assumptions, appears as an Appendix. Here we should observe that if the parameter 
B ^ 00 with A fixed then the distributional solution 

G{x) = Asgn(x), g{x) = 2AS{x) 

results, in which case for A = 1 (jl.ip becomes the Riemann shock equation 

rrit + 2{b + l)mmx = 0, 

while when B ^ with A ~ 1/5 the linear solution G{x) = Ax arises. The latter gives 
only the trivial linear PDE 

rrit + Km^ = 0, K = A rn{y) dy. 



3 Restricting the bracket to the pulsons 

To calculate the corresponding bracket for the restriction to the finite-dimensional pulson 
submanifold, we substitute m = J2Pj^i^ " Qj) i^^o both sides of ()2.2|) . which gives 

XI {{Pj^Pk}^i^ - <ij)^{y - Ik) - te, qk}Vk8{x - qj)S'{y - Qk) 

3,k 



-{Pj, qk}PjS'{x - qj)6{y - qk) + {qj, qk}PjPkS'{x - qj)S'{y - qk) 
^PjPk{G{x - y)6'{x - qj)S'{y - qk) + hG'{x - y){6{x - qj)S'{y - qk) - S'{x - qj)S{y - qk)) 



j,k 

-h^G"{x-y)5{x-q,)5{y-qk)). 
We fix e with < e < min j-Lk \qj — qk\ and then introduce the indicator functions 

Ij{x) = 1, X e [qj - e,qj + e], 
= 0, otherwise, 

and 

Jj{x) = {x- qj)Ij{x). 

Integrating the Poisson bracket relation against suitable products of indicator functions 
centred around x = qj, y = qk, such 

)Ik{y) and Ij{x)Jk{y) etc. we find 
{Pj^Pk} = -{b - l)'^G"{qj - qk)PjPk, {qj^Pk} = ip- l)pkG\qj - qk), 



{qj,qk} = G{qj - qk). (3.1) 

The brackets (jH.lj) generate the pulson 6-equations by taking the restriction of H to the finite 
submanifold i.e. 



1 f°° 1 



^Pk- 

k 

Thus we find that the pulson dynamical system p.l4p becomes 

(with G' = g). But p.ip are only Poisson brackets if the Jacobi identity is satisfied. For the 
Jacobi identity of the finite-dimensional brackets we require both 

{{qj^ qk},Pi} + cyclic = = {{pj,pk}, qi} + cyclic. (3.2) 

for all j,k,l, and similarly 

{{qj,qk},qi} + cyclic = O = {{pj,pk},pi} + cyclic. (3.3) 

The necessary conditions resulting from (j3.2j) and 13. 3|) are the functional equations (jl.lSp . 
()2.4j) . ()2.5j) . ()2.(i|) and their cyclic permutations. Therefore if G satisfies the functional 
equation then this guarantees that ()3.H1 is a finite-dimensional Poisson bracket. 



4 Conclusions 



The functional equation provides an elegant way to pick out the peakon equations 

()1.6|) from the whole class of IPDEs (ll.lj) . However, it is interesting that the compacton 
type solution G{x) = x{l — \x\/2)H{l — x)H{l + x) does not satisfy the equation. The 
compactons are solitons with compact support which appear in the Hunter-Saxton [HE] or 
Vakhnenko equations jHZllSH]- In that case G is formally a Green's function for d~^, and the 
operator 1)2.11) was formally identified as satisfying the Jacobi identity in jSHI- However, in 
[21] we noted that the Vakhnenko equation is the short wave limit of the Degasperis-Procesi 
equation, and so there may be problems with defining the boundary conditions correctly. 
Acknowledgements. We are grateful to Martin Staley for providing the Figures. Both 
authors would like to thank Toni Degasperis for introducing us to the Degasperis-Procesi 
equation. We are also grateful for the hospitality of the Isaac Newton Institute, where we 
first began working on fe = 3 in 2001. AH thanks Harry Braden and John Byatt-Smith for 
providing the general solution of the functional equation, and is grateful to Audrey Leznov 
for suggesting the simplest approach when G is odd. DH and AH also thank Simonetta 
Abenda and Tamara Grava for inviting us to the meeting Analytic and Geometric theory of 
the Camassa-Holm equation and Integrable Systems in Bologna, in September 2004, where 
we decided to finish writing up these results. DH thanks AH for taking the lead in finishing 
our project and finally transforming it from a long string of email exchanges into a finished 
product. We are both grateful to Francesco Calogero for many kindnesses and inspiring 
discussions over the years of our association together. 



Appendix: Solution of the Functional Equation by H. W. 
Braden and J.G. Byatt-Smith 

Here we shall solve the functional equation 

G'{x) [G{y) + G{z)] + G'{y) [G{z) + G{x)] + G'{z) [G{x) + G{y)] = 0, (4.1) 

where x + y + z = Q and G{x) is differentiable with continuous first derivative. We prove the 
following theorem. 

Theorem 4.1. Suppose G{x) is a solution of equation I{4-1\ ) that is differentiable with con- 
tinuous first derivative. If G{x) possesses (possibly different) power series expansions in the 
intervals (— e, 0) and (0, e) (for e> then, up to the invariance of the functional equation 

G{x) aG{Ax), (4.2) 

G{x) is one of 

(z) G{x) = - fe^" + e-^" + ll , (4.3) 
c 

which for A = gives the constant solution. 

71 

(ii) G'(x) = sin(x) sin(x ), (4.4) 

3 

which has the scaling limit G{x) = ax. 



(Hi) G{x) = sgn(x) [l - e 

We begin with some preliminary observations. 



(4.5) 



1. Equation ()4.1|) is satisfied for G a constant. Hencefortli we assume tliat G does not 
vanisli identically. 

2. G(0) G"(0) = 0. This is seen by setting = x = y = zm K^ . 

3. Suppose that both G{0) = G"(0) = 0. Then with x = and z = —y we find that 

G'{y)G{-y)+G'{-y)G{y) = 0. 

Thus G{—y)/G{y) = A, a constant. Then A = ±1, and therefore G is either even or 
odd. 

4. With X = y we obtain 

G{-2x) G'{x) + G'{x) G{x) + G'{-2x) G{x) = 0, (4.6) 

and so 



d_ 

dx 



G{-2x) 



0. 



G^{x) G{x) 
Therefore, on an interval / for which G{x) 7^ we have 

G{-2x) + 2G{x) = ci G\x) (4.7) 

It is possible for the constant cj to vary from interval to interval. This means that 
the function can have different power series expansions in different intervals. There 
are solutions possessing different power series on different intervals. For example ()4.5|1 
satisfies ()4.7j) with c = 1 for x > and c = — 1 for x < 0. (We will see how this arises 
below.) 

Thus solutions of equation (j4.7|) contain those of (j4.ip . 

5. From ()4.H) we deduce that if G' is differentiable at the point y and G{y)+G{—x—y) 7^ 0, 
then G' is also differentiable at the point x. Thus if G' is differentiable for < |xo— y| < 
e, but not at Xq we have 

G{y) + G{-x, - 2/) = 0. 
Using the continuity of G we have 

G(xo) + G(-2xo) = 0. 

We deduce from ()4.6|) that at such a point 

G'(xo)G'(-2xo) = 0. 



6. We record 

-G\-2x) = ciG{x) G'{x) - G'{x) 

2G"{-2x) = cj [G'{xf + G{x) G"{x)\ - G"{x) 
-AG"\-2x) = d [3G"(x) G"{x) + G{x) G"\x)\ - G"'{x) 
8G(™)(-2x) = C7 [4G'(x)G'"(x) + 3G"(x)2 + G(x)G(™)(x)] - G^'''\x) 



Proof of the Theorem 



The strategy of our proof is the following. We have observed that the solutions of equation 
()4.7j) contain those of ()4.H1 . We will first find the solutions of equation ()4.7|1 satisfying the 
conditions of the theorem, and then verify that these in fact solve ()4.1|) . 

Observe from ()4.7|) that if G{x) is defined in some neighbourhood of x = it may be 
extended to the whole of the real line. We will assume that G{x) has a power series expansion 
for < |a;| < e, possibly different for x positive and negative. To distinguish between the 
possibly different power series for G set 



g{x) = G{x), if X > 0, 
f(x) = G{-x), if X < 0. 



(4.9) 



Then from dHZj) 



f{2x) = c+g{xY-2g{x), (4.10) 

(7(2x) = c_/(x)2-2/(x). (4.11) 

In particular, this means that 

/(4x) = c+ [c.f{xy-2f{xf -2 [c_/(x)^-2/(x)]. (4.12) 

Our assumptions on the continuity of G and G' mean that /(O) = g{0) and /'(O) = g'{0). 
We also note from ()4.7|) that 

3/(0) = c_/(0)2, 3g{0) = c^g{0f. (4.13) 

Now we two distinct follows: 

Case (1) If /(O) = ^(0) ^ 0, then from we deduce that c_ = 3//(0) = 3/^(0) = c+. 

Therefore f{x) and g{x) satisfy the same recursion with the same initial conditions, whence 
f{x) = g{x). Together with say ()4.10j) this yields a recursion for the power series coefficients 
of the function (see for example ()4.8p ) which may be solved to give 

G(x) = - [e^^ + e"^^ + ll . (4.14) 
c 

With A = this gives f{x) = 3/c, the constant solution noted earlier. 
Case (2) If /(O) = g{0) = then there are two possibilities: 

(a) If /'(O) = g'{0) = then from our earlier observation we have that f{x) = ±g{x). 
Substituting the power series f{x) = C2X^ + . . . into ()4.10j) and employing this observation 
shows f{x) = 0, and consequently 

G{x) = 0. 

(b) Thus we may suppose /'(O) = g'{0) ^ 0. Now from ()4.8p we find that in this case 

3r(0) = (2c+-c_)/W, 
3g"{0) = {2c^-c^)g'{0f. 

1. If the second derivatives agree then c+ = c_ and again f{x) = g{x). Solving the 
recursion yields 

G{x) = sm(x) sm(x - ^ ) = — ^ Y ^ Y ~ ~6 ^ " ' ^ 

This satisfies G(— 2x) + 2G{x) = 4G^(x); that is c = 4. A scaling limit of this solution 
is 

G{x) = ax, 

which corresponds to c = 0. (This is obtained by G{x) — > aG{x), under which c — > -.) 



2. The final possibility corresponds to / and g having different second derivatives at 
X = 0. This corresponds to G' not being different iable at this point. Our earlier 
observation then yields G{ii) + G{—y) = 0. That is G is odd, whence f{x) = —g{x) 
and consequently c+ = — c_. Solving the recursion yields 

= sgn(x) [1 - e-l^'l] . (4.16) 

At this stage we have established the Theorem. 
Remark. In jS] an alternative method is described for investigating functional equations 
like (jOJ). It is instructive to compare this with the approach above. Upon setting y = x — t 
and z = t — 2x we expand ()4.H1 in t: this leads to an infinite set of differential equations 
that G must satisfy. The term independent of t yields ()4.(jj) and the equation from the order 
t term is satisfied identically as a consequence of ()4.7j) . If we substitute ()4.7|1 . which involves 
the arbitrary constant c, into the order t"^ and equations resulting from the expansion we 
can eliminate the constant c between them. This yields the following simple result 



= g{x)g"{x) [g^^^^g'{x)-g"{x)g"'' 



X 



from which we see that either g vanishes, is a linear function, or satisfies the fourth order 
differential equation given by the last term. The general solution of this latter equation is 

= ae^" + 6e-^" + ci. (4.17) 

As we have obtained a necessary condition for a function to satisfy ()4.ip we may now sub- 
stitute ()4.17|) into either the original equation or ()4.7p . placing restrictions on the constants. 
This leads to the solutions above. 

For example, we may express the functions vanishing at x = by 



_ (1 - ujx)) iujx) - P) 
uix){l + (3)c. 

where u{x) = and c+/c_ = -(1 + /5)V(1 + /3^). Then from 

^ {l-u{xr) {u{xY-(3') ^ {l-u{Ax)){u{Ax)-P) 



U(X 



l + /32)(l+/5)2c_ u{Ax){l + (3)c 



Solving this leads to /5 = and the solution ()4.5|) . ot to j3 = u, u"^, where u is a. nontrivial 
cube-root of unity and the solution (j4.4|) . (Here we have used that u{x) = e^. However, 
when j3 = 0, u, or and u{0) = 1 we may in fact deduce that u{x) = e^.) 
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